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Surface magnetic phase transition of the double-exchange ferromagnet: 

Schwinger-boson mean-field study 

Satoshi Okamotqj 
Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 

The surface magnetic phase transition of a double-exchange model for metallic manganites is 
studied using a Schwinger-boson mean-field method. About three unit-cells wide surface layers are 
identified. The magnetic moment in these layers decreases more rapidly than that in the bulk when 
the temperature is increased. This behavior is consistent with experimental observations. We also 
discuss the implication of this behavior on the tunneling magnetoresistance effect using manganites 
and possible improvement of the magnetoresistance effect near the bulk Curie temperature. 



INTRODUCTION 

Heterostructures involving transition-metal oxides 
have been attracting interest [l|, |3, y, Ij, Isj because they 
provide a playground to explore new and useful func- 
tionalities that are not realized in the bulk. In addi- 
tion, these are supposed to become fundamental building 
blocks of electronic devices utilizing a variety of proper- 
ties of transition- metal oxides [6|]. Among these oxides, 
perovskite manganites are promising candidates for spin- 
tronic devices due to their high spin polarization, high 
ferromagnetic Curie temperature (Tc), and the colossal 
magnetoresistance effect. 

One potential application of perovskite manganites is 
as a tunneling magnetoresistance (TMR) junction [Tl, |8|. 
The TMR junction consists of two ferromagnetic metal- 
lic leads separated by an insulating barrier. The con- 
ductance across the barrier can be changed by changing 
the relative orientation of magnetic moments of the two 
leads. Since high polarization in the ferromagnetic metal- 
lic leads is required to obtain the large magnetoresis- 
tace ratio (MR) , highly-polarized perovskite manganites 
could serve as ideal ferromagnetic leads. Furthermore, 
electronic devices typically function at room tempera- 
ture, therefore high Tc materials are particularly favor- 
able, such as Lai_a;Sra;Mn03 with 0.2 < x < 0.5 where 
Tc reaches 350 k[|. 

Several attempts have been made to fabricate per- 
ovskite TMR junctions 0, [HI, [H. A very large MR 
was measured at low temperature, consistent with half- 
metallicity. However, the MR decreases rapidly and dis- 
appears well below Tc [I4] ■ On the basis of spin-resolved 
photoemission spectroscopy, it was suggested that the 
rapid decrease of MR is due to the stronger temperature 
dependence of the spin polarization at interfaces than in 
the bulk [13]. 

Surface magnetism has been theoretically studied 
within a classical Heisenberg model using the numerical 
Monte Carlo (MC) technique. The surface polarization 
was shown to decrease more rapidly than that in the 
bulk [ij, [l5[ . More recently, interfacial phase transition 
of the double-exchange (DE) model for manganites was 
studied by the dynamical-mcan-ficld method [16] and the 



MC method [17[. Since the dynamical-mean- field theory 
(DMFT) neglects spatial correlations [18], it is expected 
to become less accurate in low dimensional systems, and 
therefore at surfaces and interfaces. The MC requires a 
very large system to investigate surface or interface phase 
transitions. In fact, one-dimensional systems were used 
in 17]. Therefore, the difference between the bulk mag- 
netism and the interface magnetism remains unresolved. 

In this paper, we investigate the surface magnetic 
phase transition of the DE model by using an alterna- 
tive technique, the Schwinger-boson mean-field (SBMF) 
method. We focus on metallic manganites possessing a 
relatively high Tc such as Lai_a;Sra;Mn03 with a dop- 
ing concentration of a: ~ 0.3. The SBMF method has 
had success in describing the behavior of the quantum 
Heisenberg model in low dimensions [il9] satisfying the 
Mermin- Wagner theorem [20|. The SBMF method was 
also applied to the bulk DE model [2l|, [22]. Since the 
SBMF method correctly describes low-dimensional spin 
systems, it is also expected to provide a suitable descrip- 
tion of surface and interface magnetic behavior. 

The paper is organized as follows. In section 2, the the- 
oretical model and the formalism arc outlined. In section 
3, we present the numerical results, and in section 4, we 
discuss the implication on the TMR effect and summa- 
rize. 



MODEL AND FORMALISM 

Let us start by setting up our theoretical model. Per- 
ovskite manganites are known to exhibit a variety phe- 
nomena including the colossal magnetoresistance effect, 
charge /orbital orderings and nanoscale inhomogeneity. 
In order to understand all these phenomena, many ef- 
fects must be considered. The kinetics of Cg electrons 
Hk and the Hund coupling (between the Cg electrons 
and t2g spins) Hh are the key ingredients in generating 
the DE interaction, resulting in metallic ferromagnetic 
states. In addition, there are electron-electron interac- 
tions He-e, orbital degrees of freedom, electron-lattice 
interactions He-i (Jahn- Teller coupling is also included 
here), and chemical inhomogeneity. Including all the ef- 



fects is certainly necessary to investigate the complicated 
multiphase behavior of manganites. On the other hand, 
our main interest is the surface magnetic behavior of high 
Tq manganites. Although the quantitative agreement 
between theory and experiment remains incomplete, the 
characteristic behavior of high Tc manganites can be un- 
derstood based on the simple DE model 23|, |2J]. Fur- 



thermore at this moment, detailed information on the 
surface structure is not available, although, the recent 
work suggests that the surface layer of Lai_a;Srj;Mn03 
is a (La,Sr)0 plane [25|. Therefore, we focus on the 
effect which highlights the difference between the bulk 
and the surface most: smaller coordination at the sur- 
face, by terminating the perfect cubic lattice at the [001] 
plane. For simplicity, we consider a single-orbital DE 
model, which reproduces many of the properties of metal- 
lic Lai_a;Sra;Mn03 |23j, with infinite Hund coupling. We 
will also present brief discussion on the orbital polariza- 
tion and disorder at the surface in terms of the transfer 
anisotropy. Since the theoretical model is rather simple, 
our discussion will be done on the qualitative level. 



The Hamiltonian for the DE model is written as 



H 



-tEi 



J 



jafifj + H.c.) . Here, t is the nearest- 



neighbor transfer, Sj^. the spinor boson (Schwinger boson) 
with the local constraint ^^ Si^^ia = 1 representing the 
rotation of spin space, and fi the spinless fermion rep- 
resenting an electron whose spin is parallel to the local 
moment Si = S^^g sl^SapSip. The Lagrangian for this 
system becomes 



i i<7 



E^'(-) Ei-^-i'-i ' 



(1) 



where r is the imaginary time, /i is the chemical potential 
for spinless fermions, 25 originates from the Berry phase 
of a localized spin Si, and the last term represents the lo- 
cal constraint ^^ Isjo-P = 1 with the Lagrange multiplier 

Mt). 

At this stage, we introduce a mean-field approxima- 
tion: n, = (7 J,), Xtj = (J^fj)^ and x{j = J2ai^iaSja), 

and relax the local constraint to the global one by ne- 
glecting the T dependence of A. By rescaling the spinor 

Sia, we obtain the mean- 



boson Sicr as y/zS + ni Sio 

field Lagrangians for fermions and bosons 



Lf = 22^i^^^~^^)fi^ 




H.c. 



(2) 



and 



L, 



= E **-^-*- " 9 E 



(y) 



Xij^itT^jcr 



H.C. 



E^T. 



2S 



(3) 



respectively, with 5**°* = S + ni/2. The order parameter 
X'^ is defined in the same way as before but now using the 
rescaled bosons. The local density of spinor bosons is now 
25*°'. In this SBMF, xh represents the nearest-neighbor 
ferromagnetic correlation, and the ordered moment M is 
represented by the Bose condensation of spinors Nq as 
M = iVo/2. 

We solve the self-consistency equations numerically for 
an N layer system with the open-boundary condition in 
the z direction and the periodic-boundary condition in 
the xy plane. Thus, x^ and X" are dependent on the 
layer coordinate z and the interplane bond, and A on z. 
In this mean-field theory, an additional phase transition 
appears above Tc associated with the order parameters 



X 



f,s 



Since this phase transition is an artifact of the 
present decoupling scheme, we focus on the temperature 
range below Tc. The most difficult part lies in fixing 
{Xi} so that the constraint J2cr ki^P = 25j*°* is satisfied 
at each layer. Note, J2cr kio- P at each layer depends on all 
As, and the simple bisection method does not work. Here, 
we apply the conjugate gradient algorithm and minimize 
the function A({AJ) = E/=i|E. \sia\^ ~ 25*°* |. 

As an example, we show the results of Xi and ^^ |s;ctP 
in figure [1] and order parameters x" and x^ in figure [21 
The magnetization profile corresponding to this choice 
of parameters is shown in figure [31 The Lagrange mul- 
tiplier A depends on the layer index z and temperature 
T, while the boson density remains unchanged, indicat- 
ing the accuracy of the conjugate gradient algorithm and 
applicability of the present Schwinger-boson method for 
spatially inhomogeneous systems. Typical error in the 
constraint was found to be less than 0.5 % at each layer, 
and the error in x^ and X'^ is much smaller. 

Another point to note is that our system is two di- 
mensional. Therefore, strictly speaking, there is no Bose 
condensation (no magnetic ordering) at finite tempera- 
ture unless there exists spin anisotropy. In this work, we 
discretize the momentum space and consider the lowest 
energy subband at {kx,ky) = (0,0) with the wave func- 
tion Tpo — J2a ^1=1 ^i^i<y as the Bose condensation. This 
corresponds to introducing a low-energy cutoff for spinor 
excitations representing coupling with the bulk region. In 
the following, we mainly take •\/2128 x \/2128 k points 
in the first Brillouin zone, so the typical cutoff energy 
is x^ti^k)"^ « 3 X 10~^t. The magnetic transition tem- 
peratures of the bulk system and the layer system with 
A^ = 20 were found to agree within '--^ 5 %. Thus, we 
believe that this method for a very thick system provides 
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FIG. 1: Layer dependent Lagrange multiplier A for the 20- 
layer system with a uniform transfer intensity, localized spin 
S = 3/2, and average carrier density n = 0.5. Kinks at 
T ~ 0.34t indicate the spinor Bose condensation, i.e., the 
ferromagnetic transition. The right axis shows the mean bo- 
son density at each layer, showing that the numerical error 
originating from the conjugate gradient algorithm is small. 
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FIG. 2: Order parameters xlj s-nd xli for the 20-layer DE 
model with uniform transfer. Solid lines are in-plane compo- 
nents z = Zi = Zj, while broken lines are out-of-plane compo- 
nents z = Zi — 1 — Zj. The order parameter x converges to 
the bulk value more rapidly than x" ■ 



reasonable approximation of the surface behavior. The 
wave function ■00 determines the layer-dependent spinor 
condensation (the ordered magnetic moment). 

In layered systems, n and S**"* generally depend on the 
layer index. In this work, we consider an average carrier 
density n — 0.5. Particle-hole symmetry guarantees that 
Ui and 5*°* are uniform at all temperatures. This choice 
does not lose generality as long as the system is well 
approximated by the DE model. The local spin is taken 
as S" = 3/2, therefore 5*°* = L75 

It is worth mentioning the shortcomings of the present 
SBMF method here. The bulk single-orbital DE model 
has been analyzed by using the DMFT and the ferromag- 
netic Curie temperature has been computed [26|. For 
a cubic lattice, the DMFT with the classical i2g spins 
predicts Tc/t ~ 0.2 at a; = 0.5, while the SBMF gives 
Tc/t ~ 0.36 [13. Thus, Tc is about a factor 2 overesti- 
mated in the latter. This is probably because the carriers 
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FIG. 3: Layer dependent magnetization M for the 20-layer 
DE model as a function of temperature T using the Schwinger- 
boson mean-field approximation. Left panel: surface in-plane 
transfer is taken as t%y = t, middle panel t%y = 0.5t, and right 
panel tj„ = 0.25i. 



do not suffer from scattering due to the randomly dis- 
tributed (fluctuating) spins in the SBMF. Note that the 
DMFT also tends to overestimate Tc because it neglects 
spatial correlations. In the more realistic two-band DE 
model for manganites, the on-site Coulomb interaction 
becomes one of the sources to reduce Tc [28;]. In light 
of these facts, we mainly focus on the surface magnetism 
relative to the bulk. 



RESULTS 

Numerical results for the layer dependent magnetiza- 
tion of an A^ = 20 layer system are shown in the left panel 
of figure [3] as a function of temperature. The lattice con- 
stant is taken to be unity: surface layers are located at 
z — 1 and 20, and the magnetization profile is symmet- 
ric with respect to the center of the system at z = 10.5. 
Clearly, the surface magnetization decreases faster than 
the bulk magnetization with increasing temperature, but 
all magnetizations disappear at the same temperature. 
In figure [3l many lines are on top of each other, except 
for very small z. To see the layer dependence of the mag- 
netization more clearly, we plot the magnetization as a 
function of the layer coordinate z in figure [4] for vari- 
ous temperatures indicated. One can see that layers at 
4 < z < 17 show roughly the same magnetization, thus 
about three unit-cell wide surface layers is rather robust. 
Yet near Tc, the surface layer becomes thicker. Similar 
behavior is observed in the MC study for the Heisenberg 
model [13. 

Next we discuss the effect of surface condition on the 
temperature dependent magnetization. Since we are us- 
ing the simple single-orbital DE model, we simulate vari- 
ous effects by changing the transfer intensity around the 




FIG. 4: Layer dependent magnetization M for the 20- 
layer DE model at various temperatures indicated. Filled 
(open,crossed) symbols are the results with the surface in- 
plane transfer t%y = t (0.5i, 0.25i). 



surface layers. The ferromagnetic Curie temperature re- 
mains about 0.34t in all cases, indicating that A^ = 20 
is thick enough and the surface condition does not affect 
the bulk behavior. 

First, we reduce the surface intraplane transfer tj . 
This may correspond to either the surface roughness or 
the elongation of the MnOe octahedron stabilizing the 
d^^-2_j.2 orbital in the surface layers. Results are shown 
in the middle and the right panels of figures [3] and S) 
In this case, the in-plane kinetic energy of the electrons 
is reduced in the surface layers. This causes rapid sup- 
pression of the magnetization. However, coupling be- 
tween the surface layers and the bulk region induces fi- 
nite magnetization on the surface, resulting in the long 
tail of the magnetization curve. Even in this case, layers 
at 4 < z < 17 show roughly the same magnetization. 

Second, the interlayer transfer tl between the surface 
layer and its neighboring layer is reduced. This roughly 
corresponds to the contraction of the MnOe octahedron 
in the surface layers, resulting in the increase of the 
d^2_y2 orbital occupancy. Results are shown in the mid- 
dle and the right panels of figures [5] and [H In contrast 
to the reduction of intraplane transfer, strong ferromag- 
netic correlations remain in the surface layers. This pre- 
vents the rapid reduction of the surface magnetization 
at low temperature. With increasing temperature, the 
interlayer magnetic correlation is rapidly reduced as can 
be seen in the temperature dependence of the order pa- 
rameters X* and x^ in figure [T] Thus, the surface layer 
becomes more two dimensional. Eventually, surface mag- 
netism disappears below the bulk Tc, accompanying the 
disappearance of x'^ smd x^ on the bonds connecting the 
surface layer and its neighboring layer. Above this tran- 
sition, the rest of the system behaves like the one with 
'clean' surfaces located at z = 2 and 19, since the surface 
layers are decoupled. The clear surface transition might 
be an artifact of the mean-field approximation, and in 
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FIG. 5: Layer dependent magnetization M for the 20-layer 
DE model as a function of T using the Schwinger-boson mean- 
field approximation. Left panel: uniform transfer t, middle 
panel tl — 0.5t, and right panel tl = 0.25i. 



reality the tiny magnetization on the surface layers may 
survive. If the antiferromagnetic interaction between the 
local t2g spins (neglected in the present calculations) is 
strong relative to the interplane ferromagnetic correla- 
tion, the surface magnetic moment would change its rel- 
ative orientation to the bulk moment. In either case, the 
small surface magnetic moment is expected to survive up 
to the true bulk Tc- 

The surface magnetic moment of cubic manganites was 
reported in [l3|. It was shown that the surface moment 
is much smaller than that in the bulk, but disappears at 
the bulk Tc- The experimental result seems closest to 
the theoretical curve with t^y = 0.25t in figure [3l This 
indicates that the interlayer ferromagnetic coupling re- 
mains near the surface. On the other hand in the bilayer 
manganites, the surface layer does not show a ferromag- 
netic moment, while the next-to-surface layer shows al- 
most bulk- like magnetization [29|. This situation may 
correspond to the weak interlayer coupling limit. In this 
case, the magnetization in the second layer is identical 
to that in the ideal surface. Therefore, the magnetiza- 
tion is closer to the bulk value, although it is reduced 
somewhat. But in the low Tc systems, more complexity 
would exist due to various effects such as charge/orbital 
orderings absent in the present model. 



DISCUSSION AND SUMMARY 

Here, we would like to discuss how surface (interface) 
magnetization influences the TMR effect. We consider an 
ideal tunneling junction in which an interface between a 
ferromagnetic electrode and an insulating barrier is flat. 
Furthermore, the potential barrier is very high compared 
with the band width inside the electrodes. Thus, the in- 
terface layer of an electrode is equivalent to the surface 




FIG. 6: Layer dependent magnetization M for the 20- 
layer DE model at various temperatures indicated. Filled 
(open,crossed) symbols are the results with the surface out- 
of-plane transfer tl=t (0.5t, 0.25f). 
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FIG. 7: Order parameters xtj ^-iid Xij for the 20-layer DE 
model with the surface out-of-plane transfer tl = 0.25i. Solid 
lines are in-plane components z = Zi = Zj, while broken lines 
are out-of-plane components z — Zi — 1 = Zj. 



layer in the previous discussion. The important quan- 
tity characterizing a TMR junction is the MR defined 
by MR — [Gp — Gap)/Gap with Gp(^AP) the tunnehng 
conductance with parallel (antiparallel) alignment of the 
magnetization of the electrodes [30'|. When the depen- 
dence of the tunneling matrix on the Fermi velocity v is 
weak and two electrodes are identical, MR is expressed 



MR 



2P2 



l_p2' 



(4) 



where P is the spin polarization at the Fermi level in 
an electrode defined by P = (p| — pi)/{p\ + Pi) with 
Pa the density of states of electron with spin a at the 
Fermi level. Note, MR diverges when P — 1 (full polar- 
ization). Therefore, it is not surprising to have very large 
MR for manganites at low temperature, as observed ex- 
perimentally [12| . In general, the assumptions to arrive 
at equation ([4]) are not satisfied, and one has to consider 
the realistic band structure in the presence of the inter- 
faces and the dependence of the tunneling current on the 
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FIG. 8: Temperature dependence of MR computed by equa- 
tion Q using surface magnetization (black lines, MRsurfacc) 
and bulk magnetization (gray lines, MRbuik). Solid (dashed, 
dash-dotted) lines are the results with the surface in-plane 
transfer fj^ = t (0.5f, 0.25t). Inset: the ratio between 
MRsurface and MRbuik. Experimental MRaurfacc are taken 
from U2|, and experimental MRbuik is estimated using equa- 
tion (gl) with M from [|] for a; = 0.3. 



barrier height and thickness. However, for a simple free 
electron model, it was shown that the exact solution of 
MR approaches the JuUiere's model as barrier thickness 
and height increase 3l| . Having these facts in mind, fol- 
lowing discussion will be done on a qualitative level. 

Figure [5] summarizes the results for MR using equa- 
tion (|4]). Instead of the spin polarization at the Fermi 
level, we use the total polarization P ~ M/S^°^. For 
the bulk double-exchange model, this is a rather good 
approximation. It is clearly shown that a smaller polar- 
ization at a surface layer reduces MR substantially. This 
tendency becomes stronger when the surface transfer is 
reduced. Experimental results reported in [iJl shows that 
the MR becomes as large as 800 % at the lowest temper- 
ature, but disappears far below Tc- This can be under- 
stood by the interface magnetization, which dominates 
the MR properties. Since we cannot make a direct com- 
parison between the absolute values of the theoretical MR 
and the experimental one, let us consider the ratio be- 
tween the actual MR and the MR expected from the bulk 
magnetization denoted by MRsurface and MRbuik, respec- 
tively. To compute the experimental ratio, the MR data 
in figure 3(b), [I4I, is used for MRsurface, and MR^uik is 
computed using equation ^ with Af taken from |9|] for 
the carrier concentration x = 0.3. Here, M is assumed 
to be fully saturated at the lowest temperature. The 
result is presented in the inset of figure |H1 It is shown 
that the experimental ratio MRsurface /MRbuik is compa- 
rable to the theoretical one for t% — 0.5i at temperatures 
above ~ O.STc- This may suggest that the quality of the 
experimental interfaces is rather high and the complex- 
ity inherent in doped manganites, such as charge-orbital 
ordering and (chemical) phase separation, is not so im- 
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FIG. 9: Layer dependent spectral function of spinor bosons 
computed at T = O.Oli. The small imaginary number irj with 
■q = 0.01x725*°* « 4 X 10"** is introduced in the lattice 
Green's function. Black and gray lines are the results with 
the surface in-plane transfer t%y — Q.5t and t, respectively. 
For comparison, the spectral function at z = 10 with tjj, — t 
is also shown as a gray dashed line in the lowest panel. Inset: 
a magnified view in the low energy region. 



portant near Tc- Note that in 1^ nearly optimally- 
doped manganites {x — 1/3) are used. With decreasing 
temperature, the theoretical MR increases and diverges 
at T = 0, and therefore MRsm-faco/MRbuik approaches 1. 
On the other hand, the experimental MR saturates at low 
temperature and, therefore, MRgurfacc/MRbuik decreases. 
This indicates that, in the experiment, additional inter- 
actions, such as the antiferromagnetic superexchange in- 
teraction between the localized t2g spins, compete with 
the DE interaction at low temperatures. The resulting 
canted spin structure would reduce the ferromagnetic po- 
larization and suppress the divergence of MR. 

Faster suppression of the surface magnetization in 
the theoretical results is ascribed to the larger thermal 
fluctuation of spins. To see this effect more clearly, 
we computed the layer dependent spectral function 
Az{uj) of spinor bosons. Az{uj) is defined as Az{uj) = 
-^'^J-l0j^{LJ + ir]- Hs{kx,ky)}~^\zz with uj a real fre- 
quency and irj a small imaginary number. Hs is the 
mean-field Hamiltonian for the spinor bosons given by 
the second and the third terms in L^. Az{u;) for S| and 
si have the same spectral shape. Numerical results for 
A^ = 20 systems are presented in figure [9] In the bulk 
region, x'* '^ 0.167, so the full band width of the boson 
excitation is 12x''t/2S^°^ « 0.573i. As can be seen, low- 
energy fluctuations are largely enhanced at surface layers 
z = 1 compared with those in the bulk. Even in the uni- 



form transfer case (surface tl.y = t), the low-energy part 
of the surface spectral function is about twice as large as 
that of the bulk (z = 10). Therefore, thermal excitation 
of spinor bosons has a stronger effect, resulting in the 
rapid suppression of the magnetization. It is also shown 
that the spectral function starts to develop the bulk-like 
shape at z ~ 4. This explains why surface layers are 
about three-unit-cell wide. 

The surface spectral function shown in figure [5] may 
suggest a possible way to keep a large polarization at high 
temperature to improve MR; thereby suppressing the low 
energy spin fluctuations. This may be achieved by, for 
example, (1) creating uniaxial spin anisotropy or (2) us- 
ing a magnetic insulator with a relatively high ordering 
temperature as a barrier. For (1), applying a slightly 
compressive strain for manganites would work. For (2), 
a possible candidate is BiFeOa (with [111] stacking). In 
this case, an exchange-bias-type effect is also expected. 
A parent compound of the manganites, LaMnOa, with 
[GDI] stacking may not be helpful because of its low Neel 
temperature. 

The present model includes only DE interactions. 
Therefore, at the lowest temperature, all spinor bosons 
are condensed at {kx,ky) — (0,0) and the ferromagnetic 
moment is saturated. In reality, additional interactions 
may create complexity. For example, surface polarity 
and segregation [25|, |32] are expected to change the sur- 
face carrier density and reduce the ferromagnetic inter- 
action. Therefore, antiferromagnetic superexchange in- 
teractions between the localized t2g spins are expected 
to reduce the ferromagnetic polarization at low temper- 
ature. Further complexity, such as polaronic effects and 
charge-orbital ordering, would enhance the effect of the 
superexchange interactions. The present SBMF method 
is rather simple, and including such complexities to the 
fermionic part is possible. We are currently working on 
including the electron-lattice couplings with orbital de- 
generacy using a dynamical-mean-fleld-type treatment to 
study the magnetic and metal-insulator transitions at the 
surface and in thin films [33|. 

To summarize, we studied the surface magnetic be- 
havior of the double-exchange model for doped mangan- 
ites using the Schwinger-boson mean-field method. Low- 
energy spin fluctuations are enhanced at the surface and 
the magnetic moment is suppressed more rapidly than in 
the bulk. We further considered an ideal tunneling mag- 
netoresistance junction consisting of two manganite leads 
and an insulating barrier. The magnetoresistance ratio of 
such a junction is determined by the surface polarization. 
Therefore, it decreases much faster than the bulk mag- 
netization when the temperature is increased. A possible 
improvement of the MR is expected from suppressing the 
low-energy spin fluctuations. 
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